Electronic states and superconductivity in multilayer high-T c cuprates 
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We study electronic states of multilayer cuprates in the normal phases as functions of the num- 
ber of CuC>2 planes and the doping rate. The resonating valence bond wave function and the 
Gutzwiller approximation are used for a two-dimensional multilayer t-t' -t" -J model. We calculate 
the electron-removal spectral functions at (7T,0) in the CuC>2 plane next to the surface to understand 
the angle-resolved photoemission spectroscopy (ARPES) spectra. We find that the trilayer spectrum 
is narrower than the bilayer spectrum but is wider than the monolayer spectrum. In the tri- and 
tetralayer systems, the outer Cu02 plane has different superconducting amplitude from the inner 
CUO2 plane, while each layer in the bilayer systems has same amplitude. The recent ARPES and 
NMR experiments are discussed in the light of the present theory. 
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I. INTRODUCTION 

In studies of high-T c superconductivity, the number 
of Cu0 2 planes in a unit cell, n, is one of the key factors. 
It is indicated that the maximum T c monotonically in- 
creases with n and saturates at either n—3 or 4. 1 ' 2 The 
bilayer compounds have crystallographically equivalent 
Cu02 planes, while the tri- and tetralayer compounds 
have inequivalent planes, i.e., the outer pyramidal coor- 
dinated planes (OP) and the inner square coordinated 
planes (IP). By calculating the Madelung energy for Bi 
and Tl compounds in the point-charge model, it was 
shown that the hole density of the IP, A^(IP), is lower 
than that of the OP, AT^(OP), since the IP is sandwiched 
by Ca sheets with positive charge. 3 ' 4 It was also argued 
that a hole density difference AN h = N h (OP) - N h (IP) 
becomes large by increasing n and controls the supercon- 
ducting phases. 5-7 

Such hole distribution among the Cu02 layers 
was observed above T c by the 17 and 63 Cu nu- 
clear magnetic resonance (NMR) measurement for Pb- 
doped Bi 2 Sr 2 Ca2Cu30io +y (Bi2223), 8 < 9 and for Hg 
Sr 2 Ca 2 Cu 3 Oio+y (Hgl223). 10,11 These studies showed 
that the IP has smaller carrier density and larger antifer- 
romagnetic (AF) correlation than the OP. Kotegawa et 
al. systematically estimated the hole density in each layer 
by the 63 Cu Knight shift in HgBa 2 Ca n _iCu n 2n + 2 +y 
and CuBa 2 Ca n _iCu n 2n +4_y. 12 By use of an empirical 
relation between the Knight shift at room temperature 
and the hole density, which was deduced from the nuclear 
quadruple resonance (NQR) frequency, 13 they found that 
the hole density difference among the layers increases as 
either a total carrier content Stotai or n increases. 

As regards an interlayer hopping, the bilayer split- 
ting was shown by band calculations 14-16 and clearly ob- 
served by the angle-resolved photoemission spectroscopy 
(ARPES) in Bi 2 Sr 2 CaCu 2 8+y (Bi2212) above T c . 17-19 
However, a trilayer splitting in Bi2223 has not been 
observed. 20,21 To understand the electronic states in the 
multilayer cuprates, it is necessary to clarify the rea- 



son why the trilayer splitting is not clearly observed, al- 
though the bilayer splitting is distinguishable. 

In the present paper, the electron-removal spectral 
function is calculated in a two-dimensional multilayer t- 
i -i"-,/ model at zero temperature by using the resonating 
valence bond (RVB) wave function and the Gutzwiller 
approximation. We study the n dependence of the 
single-particle excitations including the interlayer hop- 
ping which connects the chemically different Cu0 2 planes 
for n > 3. Since the ARPES measurement is sensitive to 
a surface of sample due to a short escape depth of outgo- 
ing electron, the spectra in the tri- and tetralayer com- 
pounds should be dominated by the OP's contribution. 
Moreover, the spectra depend on the incident photon en- 
ergy hv and the photoemission matrix element. 22-25 The 
spectra in the bilayer systems have two peaks composed 
of the bonding (B) and anti-bonding (AB) bands, while 
the spectra in trilayer systems have three peaks. We find 
that the splitting between the two dominant peaks in 
the trilayer system is smaller than the bilayer splitting. 
Then, the trilayer spectrum is narrower than the bilayer 
spectrum but is wider than the monolayer spectrum. 

In the superconducting (SC) phase, Tokunaga et 
al. studied the temperature dependence of the 
Knight shift {K s ) and the nuclear spin-lattice relax- 
ation rate (1/Ti) of 63 Cu in the tetralayer compound, 
(Cuo.6Co.4)Bi 2 Sr 2 Ca3Cu40i 2+y , 26 and found two char- 
acteristic temperatures associated with the IP and the 
OP, T (IP) and T (OP). Both K s and l/T x in the IP 
markedly decrease below Xb(IP) = 117 K, while those in 
the OP moderately decrease below Tb(IP)=117 K, but 
markedly decrease below To(OP)=60 K. These results 
suggest that a bulk SC transition set in at 117 K but the 
temperature dependence of SC amplitude in the OP is 
different from the BCS theory. Such a non-BCS ampli- 
tude is found in a two-band system where the SC am- 
plitude in one band is different from that in the other 
band, and these bands are coupled to each other by a 
pair-transfer interaction. 27 The two kinds of Tq's are ob- 
served in the tetralayer compounds, while it is not found 
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in the trilayer compounds. 12 We calculate the SC ampli- 
tude of each layer and find that the outer CuC>2 plane 
in both the tri- and tetralayer systems has different SC 
amplitude from the inner CuC-2 plane. In our calcula- 
tions, the difference of SC amplitude in the trilayer sys- 
tems is comparable to that in the tetralayer systems. We 
consider that the To(IP) different from Tb(OP) could be 
observed also in the trilayer compounds. 

The rest of this paper is organized as follows. In Sec. 
II, we introduce the multilayer t-if-if'-J model that in- 
cludes the interlayer hopping and the site potential. Sec- 
tion III is devoted to show the electronic states in the nor- 
mal phases and the electron-removal spectral functions. 
The layer dependence of superconducting amplitudes is 
shown in Sec. IV. Summary and discussions are given in 
Sec. V. 



II. THEORETICAL MODEL 

We examine the electronic state in the characteristic 
block of n CuC-2 layers. The Hamiltonian of each layer 
is given by 



H intra — t s C„, „■ „C, 



a,(ij)iat,<r a,{ij) 2 „d^ 

-*" E + H.c. 



(i) 



a,(ij)i 



where c a .i M — c ai j jCT (l — n ai j_j) is the annihilation 
operator of an electron in the a layer with spin a at 
site i with the constraint of no double occupancy, and 
n at i-„ and S a ^ are the charge and the spin operators, 
respectively. The summations (ij)i st , (ij)2nd, and (ij) 3r d 
run over first, second, and third nearest-neighbor sites, 
respectively. The values of the parameters are as fol- 
lows; J = 0.14 eV, t/J = 2.5, t'/J = -0.85, and 
t"/J = 0.575. 28 ' 29 

The interlayer hopping has the dispersion relation 14 



e±(k) = \cos(k x ) - cos(fc y )] 2 , 



(2) 



which is obtained by integrating out the high-energy de- 
grees of freedom from the Hamiltonian reproducing the 
full local-density approximation (LDA) calculation. 15,16 
We adopt Eq. (2) for n > 2 and t±/J = 1.0, 14-18 and 
write the Hamiltonian of the interlayer hopping as 

H inte r= E £ -l( /c ) (Cfc,^,/^ + H.c.) , (3) 

where a and f3 are the neighboring layer indices. 

To describe the hole distribution among the layers, we 
introduce a site potential W in the IP's as 



H s ite = W E ^ l k <J c ltk ,a = W E n 7 , 



(4) 



7— in,/c,<T 



where 7 indicates the IP index. The value of W should 
be negative. For simplicity, W is taken to be independent 
of the hole densities. 

The total Hamiltonian is 



H = Hi, 



~t" Hinter ~t" H. 



site • 



(5) 



Since this Hamiltonian includes the constraint of no dou- 
ble occupancy, we adopt the Gutzwillcr approximation 
and statistically average the constraint before variational 
calculation. Details of calculations are summarized in the 
Appendix. In the normal phases, the order parameter, 



is introduced and referred to as "uniform RVB" (URVB). 
In addition to this, in the SC phase the order parameter, 



B n 



T jy ^ ] (Cq,t,|Cq ; i-)-T,J. C Ct ,», J. Cq ,i+T, ] ) ) (7) 



is included and referred to as "singlet RVB" (SRVB). 

Then, we obtain the following mean-field Hamiltonian, 
which is numerically solved: 

Hrvb - H N 

+ E 

- E ~J S* )T cos(fe T )c 



k,r—x,y 



Xct.i 



+B a . T cos(k T )c a k ^c a _ k 
-^2t[cos(k x ) + cos(k y )} + At' cos (k x ) cos(k y ) 

+2i" [cos(2k x ) + cos(2fcj,)] } - fj, 

3 » 

-~;J[Xa,x cos(fc 

x) "b Xa,y COs(fcy)] , 

= 4EE cos ( fc ^ c lfc,- c «.^>' 



(8) 



k a 



~ E 2cos(fc r )(c QiM c Q -k,i), 



(9) 
(10) 

(11) 



where /j, is the chemical potential and N is the total 
electron number. The values of hopping integrals and 
the magnetic coupling constant are renormalized by the 
Gutzwiller factors as 

t = g t ■ t, t' = g t ■ t', t" = g t ■ t" , t± = g t ■ t±, 

J = 9J- J, (12) 
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FIG. 1. W dependence of AN h in the URVB state for the 
trilayer and the tetralayer with S = 0.25 and 0.3. The values 
of the parameters are J = 0.14 eV, t/J = 2.5, t'/J = —0.85, 
t" /J = 0.575, and t±/J = 1.0. The averaged doping rate 5 is 
defined as 8 — 5 to tai/n, in which S to tai is a total hole density 
in a unit cell. 



and 



g t = 26/(1 + 6), 5J = 4/(1 + 5) 2 



(13) 



The even parity is imposed on Xa,r- The k summations 
are carried out on 160x160 points in the first Brillouin 
zone. 



III. ELECTRONIC STATES IN NORMAL PHASES 

We study single-particle excitations to clarify the elec- 
tronic states in the normal phases. The characteristics of 
multilayer systems are manifest at (tt,0) in the dispersion 
relations, since the interlayer splitting is largest at this 
point according to Eq. (2). The electron- removal spectral 
functions at (n,0) are shown in this section. The aver- 
aged doping rate 6 is defined as <5 = Stotai /n, in which 
6totai is a total hole density in the characteristic block of 
n CuC>2 layers. We choose 6 = 0.25 or 0.3, below. These 
values of S are useful to study the n dependence of the 
single-particle excitations, since the ARPES spectra at 
(7r,0) show sharp quasiparticle peaks for a large doping 
rate. On the other hand, the ARPES spectra become 
broad and dull by decreasing the doping rate. ' 31 

Before proceeding to the study of ARPES spectrum, it 
is necessary to determine the value of W that is still a free 
parameter. We calculate ANh as a function of W in the 
URVB state of tri- and tetralayer systems with 5 =0.25 
and 0.3. The results are shown in Fig. 1. We determine a 
value of W by comparing Fig. 1 to the results obtained by 
Kotegawa et al. 12 in NMR measurement. Their results of 
ANh in the tri- and tetralayer samples are fitted by the 
following equations: 

trilayer, AN h = -0.130 + 0.759 6; (14) 
tetralayer, AN h = -0.150 + 0.934 S. (15) 



By using Eqs. (14) and (15), the values of W that repro- 
duce the NMR results are determined as 



trilayer, 

r 6 = 0.25 
I 5 = 0.30 
tetralayer, 

f 5 = 0.25 
\6 = 0.30 



AN h 
AN h 

AN h 
AN h 



0.06 
0.10 

0.08 
0.13 



W/J 
W/J 

W/J 
W/J 



-0.25, 
-0.35; 

-0.30, 
-0.45. 



(16) 
(17) 



The dispersion relations of tri- and tetralayer systems 
are plotted together with the bilayer system (W = 0) 
in Fig. 2 for 5 = 0.25. Since the interlayer hopping 
is given by Eq. (2), the dispersion relations in the bi- 
layer system are degenerate in the (0,0)-(7r,7r) direction 
and the splitting becomes maximum at (7T,0). The band 
with lower-binding energy indicated by the solid line in 
Fig. 2(a) is the AB band and the other one is the B band. 
The bilayer splitting at (n,0) for 6 = 0.25 is 0.8J, which 
is given by 2t±g t . Most of experiments 17-19,22,23 and 
band calculations 15 ' 16 show two holelike Fermi surfaces 
consistent with our dispersion relation. On the other 
hand, there are experimental 32 and theoretical studies 33 
suggesting an clcctronlike Fermi surface, where the AB 
band crosses the Fermi level between (0,0) and (7T,0). 
Such a dispersion relation would be achieved by small 
changes of parameters. In the trilayer system, there ex- 
ist three bands. Among them, one of the bands indi- 
cated by the dotted line in Fig. 2(b) is separated from 
the other two in the whole Brillouin zone due to the 
site potential. The middle band shown by the solid 
line in Fig. 2(b) is composed of only the OP's degrees 
of freedom as, (|OPi) - |OP 2 ))/V2, where |OPi) and 
IOP2) indicate the wave functions in the outer planes. 
We refer to this band as the "antisymmetric (AOP) 
band" . The wave functions of remaining two bands are 
/(W)(|OPi) + \OP 2 ))/V2±g(W)\IP), where |IP) indi- 
cates the IP's wave function. The coefficients f(W) and 
g(W) are functions of W. The positive (negative) sign 
corresponding to the dotted (broken) line in Fig. 2(b) is 
referred to as the "IP band" ("symmetric (SOP) band"). 
In the tetralayer system, OP's and IP's combine to give 
symmetric and antisymmetric combinations as |OP±) = 
flOPr) ± |OP 2 ))/V2 and |IP±) = (|IPi> ± |IP 2 ))/x/2. 
The four bands shown in Fig. 2(c) have the follow- 
ing combinations at (jr, 0) from the top to the bottom, 



I a) = /i|IP_) + ff i|OP_), 



/ 2 |IP+) - 52 |OP+), 



|7> = /s|IP-) -5 ? |OP_), and \5) = / 4 |IP+) + 54 |OP + ), 
where fi and gi [i = 1 ~ 4) are functions of W. 

In the following, we assume that the ARPES experi- 
ments in Bi compounds mainly detect the single-particle 
excitation in the Cu0 2 plane next to the BiO layer. The 
ARPES measurement is sensitive to a surface of sample 
due to a short escape depth of outgoing electron, which 
is about 10 A or less depending on its kinetic energy in 
the range of 20—200 eV. 34 The cleaved surface is consid- 
ered to be the Bi-0 layer, from which the distance to the 
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FIG. 2. Dispersion relations of (a) bilayer, (b) trilayer, 
and (c) tetralayer systems for 8 — 0.25 in the URVB state. 
The values of the parameters are, J = 0.14 eV, t/J = 2.5, 
t'/J = -0.85, t"/J = 0.575, and t±/J = 1.0. 
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FIG. 3. The n dependence of the electron-removal spectral 
function in the normal phase for J = 0.14 eV, t/J — 2.5, 
t'/J = -0.85, t"/J = 0.575, and t ± /J = 1.0. The compo- 
nents in the Cu02 plane nearest to the surface are plotted (a) 
on the n-uj plane for 8 =0.25, and (b) as functions of uj for 
8 =0.30 for the bi- and trilayer systems. The thin solid (bro- 
ken) curve is obtained by performing a Lorentzian broadening 
with a width of 0.35J for the 8 functions in the bi- (tri-) layer 
system. In (a) and (b), the values of W follow Eqs. (16) and 
(17). 



IP of trilayer compounds is also about 10 A. 1,2 Therefore 
the ARPES spectra in the tri- and tetralayer compounds 
should be dominated by the OP's contribution. More- 
over, the spectra depend on the incident photon energy 
hv and the photoemission matrix element. 22-25 In the 
case of overdoped bilayer samples, an antibonding (AB) 
band has a sharp peak at (ir,0) for hv ~ 22 eV. Although 
a bonding (B) band is unclear for hv ~ 22 eV, it becomes 
clear for hv - 32 eV (Ref. 23) and 47 eV. 22 

The n dependence of the OP contributions in the 
electron-removal spectral function A~(k,Lu), at (7r,0) is 
plotted in Fig. 3 for t±/J = 1.0. The energy is mea- 
sured from the Fermi level. In the bilayer system, the two 
quasiparticle peaks clearly separate each other for both 
doping rates. The peak with lower-binding energy is the 
AB band and the other is the B band. In the trilayer sys- 
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TABLE I. The binding energies and the ARPES spectral 
weights at k — (it, 0) in a trilayer compound for <5=0.3. Each 
contribution from the IP and the OP is listed. 





LU/J 


OP" 


IP 6 


SOP-band 


0.14 


0.31 


0.37 


AOP-band 


0.64 


0.50 


0.00 


IP-band 


1.49 


0.19 


0.63 


a A C11O2 layer 


nearest to the surface. 






b A CuC>2 layer 


next nearest to the surface. 




TABLE II. The binding energies 


and the ARPES spectral 


weights at k - 


= (n, 0) in a tetralayer compound for S 


=0.3. 


Each contribution from the IP and the OP is listed. 






lo/J 


OP a 


IP 6 


a-band 


0.08 


0.25 


0.25 


/3-band 


0.35 


0.43 


0.07 


7-band 


1.00 


0.25 


0.25 


5-band 


1.65 


0.07 


0.43 


a A C11O2 layer 


nearest to the surface. 






6 A CuC>2 layer 


next nearest to the surface. 





tem, the OP's combine to give symmetric and antisym- 
metric combinations, the latter of which is the AOP band 
as denned above and is located at u> / J— 0.64. Only the 
former can hybridize with the IP, and produces the SOP 
band (cj/J=0.14) and the IP band (w/J=1.49). The to- 
tal splitting in the trilayer compound (Aw/J=1.35) is 
larger than in the bilayer compound (Aw/J=0.9), but 
the SOP- AOP splitting (Aw/J=0.5) is smaller than the 
bilayer splitting. Therefore the trilayer spectrum seems 
to be narrower than the bilayer spectrum as shown in 
Fig. 3(b), in which the appropriate broadening is given. 

The contributions to the spectral weight from the OP 
and IP are listed in Table I for the three bands in the 
trilayer system together with their binding energies. The 
OP contribution distributes among the three bands as 
SOP:AOP:IP-3:5:2, while the IP contribution is ob- 
tained as SOP:AOP:IP~4:0:6. If the electron escape 
depth is longer than 10 A , the weight of the IP band in- 
creases as a result of the IP contribution to the ARPES 
spectrum. This makes the spectrum broader than the 
bilayer spectrum as easily expected from the Fig 3(b). 
Thus, the broadness might be a measure of the escape 
depth. 

Each spectral weight in the tetralayer compound is 
listed in Table. II. The dominant three peaks in the OP, 
i.e., a, (3 and 7 band, may result in a broad peak due to 
the line overlap. Including the IP contribution, it may 
become more difficult to distinguish the peaks. 



IV. SUPERCONDUCTING AMPLITUDES 

In this section, we consider the superconducting (SC) 
state, where the SRVB state is included in addition to 
the URVB state. As observed in the tetralayer com- 
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FIG. 4. The n dependence of the SC amplitude at S = 0.3. 
For n > 3, the up and down triangle indicate the amplitude 
in the OP and the IP, respectively. At n =2, both layers have 
the same amplitude. The value of parameters are, J = 0.14 
eV, t/J = 2.5, t'/J = -0.85, t"/J = 0.575, and t±/J = 1.0. 
The values of W follow Eqs. (16) and (17). 



pounds by NMR experiment, the IP and the OP have 
different Tq's, which are ascribed to the different hole 
concentration in each layer. 26 In the trilayer compounds, 
the different To's are not still found, 12 although there is a 
finite ANh , which is smaller than ANh in the tetralayer 
compounds. To understand the SC state in the multi- 
layer cuprates, we calculate the SC amplitude in each 
layer defined as \gtB a , T \, and show the n dependence. 
In the following, we assume that all Cu02 layers are in 
the SC phase, and the value of W in the SC phase is 
the same as that in the normal phase. The results are 
shown in Fig. 4. The SC gap has the d x 2_ y 2 symmetry 
in each layer for any n. The OP's in both the tri- and 
tetralayer systems have different superconducting ampli- 
tude from the IP, while each layer in the bilayer systems 
has the same amplitude. For n > 3, the difference of SC 
amplitude slightly increases with increasing n. Since the 
difference in the trilayer systems is comparable to that in 
the tetralayer systems, the different To could be observed 
in experiment. 

We note that, although the SC amplitude is larger in 
the OP than the IP as seen in Fig. 4, the relative magni- 
tude depends on the intralayer's parameters, which mod- 
ify the phase diagram of ground state in the monolayer 
system. 35 In the present calculation, moreover, the site 
potential is fixed at the value in the normal phase. It 
is also possible that the site potential depends on the 
hole densities and is determined self-consistently. If such 
mechanisms are taken into account, the hole distribution 
among the layers may changes below T c and the valance 
of SC amplitudes can be modified. These possibilities are 
studied elsewhere. 

We did not mention quasiparticle excitations in the SC 
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phase. The coherent SC peak in Bi2212 seems to split 
into two peaks. 17 There are some interpretations of the 
bilayer splitting in the SC phase. 17,33 To clarify its origin, 
more experimental studies on the multilayer cuprates are 
needed. 



V. SUMMARY 

We have studied the electronic states and the super- 
conductivity in multilayer cuprates by use of the resonat- 
ing valence bond wave function and the Gutzwiller ap- 
proximation for the two-dimensional multilayer t-t'-t"- 
J model. In this model, the one-particle hopping con- 
nects the chemically different Cu02 planes in the tri- and 
tetralayer systems. We calculated the electron-removal 
spectra in the Cu02 plane next to the BiO layer, since 
the ARPES is sensitive to the surface of sample. We find 
that the ARPES spectra in the trilayer systems at (ir,0) 
should be narrower than that in the bilayer systems, and 
wider than that in the monolayer systems. This is con- 
sistent with recent experimental data. 36 

The superconducting (SC) amplitude of each layer is 
shown as a function of n. In bilayer system, both layers 
have the same SC amplitude. In tri- and tetralayer sys- 
tems, the SC amplitude of the OP is different from that 
in the IP due to the different hole concentration. Such 
a difference could be observed as two different T also in 
the trilayer compounds. 
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APPENDIX: GUTZWILLER APPROXIMATION 
AND RVB WAVE FUNCTION IN MULTILAYER 
SYSTEMS 

The RVB-mean-field theory with the Gutzwiller ap- 
proximation showed that the d x 2_ y 2 wave supercon- 
ductivity can be stabilized in the two-dimensional t- 
J model. 37 The Gutzwiller approximation has been re- 
cently generalized to the coexistent state of the d-wave 
SC and AF order. 38 In this section, we summarize the 
process to treat the constraint within this approxima- 
tion. 

Explicitly describing the constraint, the Hamiltonian 
Eq. (5) has the following relation to the nonconstrained 
Hamiltonian: 



where Pq is the Gutzwiller's projection operator and de- 
fined as 

P G = l[(l-h n h il ). (A2) 

i 

We assume the projected RVB wave function 

\il>)=PG\MXc,T,B a , T ,n)), (A3) 

where Xu,t and _B QjT are the variational parame- 
ters relating to URVB and SRVB, respectively, and 
\ipo(Xa,r> B a>T , [/,)) is a mean-field wave function with 
URVB and SRVB orders. Using these notations, the vari- 
ational energy E var = (H) is rewritten as E v&1 = (H)q, 
where the parameters t and J in H are replaced with 

t = gt- t, t' = g t ■ t', t" = g t ■ t", t± = g t ■ t±, 

J = 9J-J, (A4) 



and 



gt = 26/(1 + 5), gj = 4/(l + 6) 2 



(A5) 



In the Gutzwiller approximation, the effects of the pro- 
jection are statistically averaged and renormalized into 
the expectation values as follows: 



( C L C ^> = 9t{c] a c ja ) , (S i -S j }=gj(S i -S j )o, (A6) 

where (• • -)o represents the expectation value in terms of 
\ipo) = \ipo(Xa,T, B aiT , n)), and (• • •) represents the nor- 
malized expectation value in \tp) = Pg\^)o\ 



- n.lA 



(0) = 



\0\ 



(V»o|P G 0P G |V>o) 



(A7) 



H = P G HP G , 



(Al) 



By considering the rcnormalization factors of coupling 
constant, the remaining work is to solve the self- 
consistent equations, in which the renormalization fac- 
tors are also determined self-consistently. In the multi- 
layer systems, the constraint is included in the interlayer 
hopping and renormalizes t±. Eder et al. numerically 
showed that the t± is renormalized by an in-plane quasi- 
particle weight. 39 Such an effect is reproduced by the 
Gutzwiller factor g t . 

Here, we summarize the eigenvalue problem in the SC 
state of multilayer systems. 40 Neglecting the details, the 
Hamiltonian generally has the following form: 

H = Ac l,k,* C P-,k,* 

= <^ M( t> + Tr A ( A8 ) 
k 

4> ] = (cL, T ' c «.-m) . (A9) 



G 



where 



^ = A, t B = B. 



(A10) 



The matrix M is hermite and satisfies the following rela- 
tions: 



a 3 criMaia 3 = —M*, 
ikft = M. 



(All) 
(A12) 



We obtain the eigenvalue and the eigenvectors by solving 

X r . 



MV n =LO n V n = UJ n 



(A13) 



By taking account of Eq. (All), the complex conjugate 
of Eq. (A13) satisfies the following equation: 



Ma!a 3 V* = -u} n a-iazV*. 



(A14) 



Therefore W n = a\a 3 V* is also the eigenvector of M. 
If V n belongs to a positive eigenvalue u> n > 0, another 
eigenvector W n belongs to a negative eigenvalue —u> n . 
By use of bj, and b n defined below: 



i, \ / v* v* 

u n \ \ I -^n.ct 1 n,a 



= E 



Y X I \ l A 



(A15) 



the Hamiltonian Eq. (A8) is diagonalized as 

I? = 2 5>„6t 6„ - + TrA (A16) 

n>0 n>0 

By using Eq. (A15), the electron operator is 

c Q , T = ]T - Y* in %) . (A17) 

n 

The order parameters and the spectral weight are calcu- 
lated as 

(o|4 iCT c^|o) = ^y; i „y„, a , (ais) 

n 

(oi4, CT 4,-ai°) = E-^« y ™.- ( A19 ) 



|(n|c QiCT |0)| 2 = \Y a 



(A20) 
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